Orientation-dependent spontaneous emission rates of a two-level quantum emitter in 

any nanophotonic environment 
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We study theoretically the spontaneous emission rate of a two-level quantum emitter in any 
nanophotonic system. We derive a general representation of the rate on the orientation of the tran- 
sition dipole by only invoking symmetry of the Green function. The rate depends quadratically on 
orientation and is determined by rates along three principal axes, which greatly simplifies visual- 
ization: Emission-rate surfaces provide insight on how preferred orientations for enhancement (or 
inhibition) depend on emission frequency and location, as shown for a mirror, a plasmonic sphere, or 
a photonic bandgap crystal. Moreover, insight is provided on novel means to "switch" the emission 
rates by actively controlling the orientation of the emitters' transition dipole. 



I. INTRODUCTION 

It is well-known that the characteristics of sponta- 
neously emitted light depend strongly on the environ- 
ment of the light source [J HSU- According to quan- 
tum electrodynamics, the emission rate of a two-level 
quantum emitter, described by Fermi's golden rule, is 
generally factorized into a part describing the sources 
intrinsic quantum properties and another part describ- 
ing the influence of the environment on the light field. 
Currently, there are many efforts to control the emission 
rate of quantum emitters by optimizing the nanoscale 
environment by, e.g., reflecting interfaces IlLSj, micro- 
cavities d, 0] , photonic crystals [1, [t| [l(| HI | , or plas- 
monic nanoantennae [HI, O EH • Control of spontaneous 
emission is notably relevant to applications, including 
single-photon sources for quantum information, minia- 
ture lasers an d lig ht-emitting diodes, and solar energy 
harvesting [IH [la [l7l| . 

The effect of the environment of a source on its emis- 
sion rate is described by the local density of optical states 
(LDOS) [l,i,[iH. The LDOS counts the number of pho- 
ton modes available for emission, and it is interpreted 
as the density of vacuum fluctuations. In many exper- 
imentally relevant cases, it is theoretically known that 
emission rates strongly differ for various orientations of 
the transition dipole moment see, e.g., [1, Ell- Thus, the 
widely pursued control of position and frequency leaves 
a large uncertainty in the emission rate To date, no 
clear picture has emerged of the general characteristics 
of the orientation dependence. It is an open question 
whether the behavior mimics the local symmetry around 
the emitter, see Fig. [TJ or whether any generic depen- 
dence exists at all. 

Therefore, we present fundamental insights in the com- 
plex dependence of the emission rates of a quantum emit- 
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FIG. 1: (color online) Drawing of a two-level quantum emit- 
ter embedded in an arbitrary nanophotonic system, here de- 
picted as a cluster of 6 scatterers. If the emission rate were 
to mimic the symmetry of the system, one would here expect 
an emission-rate surface with a 6-fold symmetry. Our analy- 
sis reveals, however, that these surfaces take on only specific 
shapes determined by the symmetry of the Green dyadic. The 
symmetry analysis allows one to conclude without any calcu- 
lation that the rate is identical for all dipole orientations in 
the plane of the 6 scatterers. 



ter on the orientation of its dipole moment. Our general, 
yet simple theoretical analysis only invokes the symme- 
try of the Green function and provides a complete classi- 
fication of the orientation-dependences that the emission 
rate can assume in any nanophotonic system. This classi- 
fication leads to an intuitive visualization that is based on 
only a few clearly defined physical parameters, as shown 
by examples of an emitter near a mirror, a plasmonic 
sphere, or in a 3D photonic bandgap crystal. From our 
analysis, we conclude that control over the orientation of 
the transition dipole moment opens novel applications: 
If one can tune the orientation of an emitter, one can 
"switch" emission from inhibited to enhanced and vice 
versa. In the field of quantum information [l9l ]. atomic 
qubits that fly by nanophotonic systems could acquire 
controllable phase shifts by tuning their orientation rel- 
ative to the principal axes. 
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II. THEORY 
A. Derivation of emission-rate surface 

The rate of spontaneous emission T of a two-level dipo- 
lar quantum emitter in the weak-coupling approximation 
is equal to [I,©, [Tl|: 

r(w, r, e d ) = A(r, u, e d ), (1) 

with u) the emission frequency, r the source's position, 
the dipole orientation, d the modulus of the matrix 
element of the transition dipole moment. N(r,uj,e d ) is 
the local density of optical states (LDOS) that equals: 

N(u, r, ej) = -^(ej • Im(G(r, r, u)) ■ e d ), (2) 

7TC 

with G(r, r, lo) the Green dyadic Q. Eq. (fl) reveals the 
well-known fact that the emission rate depends on the 
frequency and the position of the emitter. As is well 
known, Eq. @ is also applicable to emission dynamics 
inside dissipative optical media. In such media, the imag- 
inary part of the Green dyadic describes the total decay 
rate, i.e., the sum of the radiative decay rate and the rate 
of quenching induced by the environment. Hence, the re- 
sults in this paper carry over straightaway to the decay 
dynamics of dipolcs emitters in dissipative nanophotonic 
environments. 

A didactic example to illustrate the dependence of 
emission rates on frequency, position and dipole orien- 
tation is that of a source near a perfect mirror, see 
Fig. [^AY, which can be understood from image dipole 
analysis [l|, [3j ■ The rate depends strongly on the dipole 
orientation e^: at small distances a dipole parallel to 
the mirror has a vanishing emission rate, which can 
be interpreted as due to destructive interference of the 
dipole with its oppositely oriented image. In contrast, a 
dipole perpendicular to the mirror has twice the unper- 
turbed rate owing to constructive interference, as shown 
in Fig. [U^B) ■ Clearly, the symmetry of this particular ge- 
ometry implies that the parallel and perpendicular dipole 
orientations are 'principal' orientations along which the 
maximum and minimum rates are attained. At interme- 
diate orientations the rate is a weighted average of the 
two rates. 

The main result of our paper is that the rate always 
depends on orientation via a quadratic form with three 
perpendicular principal axes, as will now be proven: On 
account of reciprocity, the Green dyadic is equal to its 
transpose upon exchanging the coordinates. Hence 

Im(G(r, r', uj)) t = Im(G(r', r, u)). (3) 

Furthermore, the imaginary part of the Green dyadic is 
real. Therefore, the imaginary part of the Green dyadic 
in Eq. Q is a real and symmetric 3x3 matrix. Conse- 
quently, at each frequency to and spatial position r, the 



imaginary part of the Green dyadic can always be diago- 
nalized, and has 3 eigenvalues (gi, <?2, 33) that correspond 
to three orthogonal eigenvectors. Since the eigenvalues 
can be ordered by magnitude, we relabel the eigenvalues 
and the concomitant main axes as {v m j n , v mc( j, v max }. 
This basis corresponds to three perpendicular principal 
dipole orientations that vary with dipole location r and 
frequency u>. In this orthonormal basis we express the 
dipole orientation unit vector e d as: 

ed = /3lV min + /3 2 V mc d + AjVmax, (4) 

where Pi are coefficients that are constrained through 
1 +02 + 0i = 1 to lie on a unit sphere, since He^H = 
1. Clearly, the coefficients are functions of the dipole 
orientation: = 0(e d ) = ej ■ v 4 . 

Using Eqs. Hid]), the emission rate T can be expressed 
in emission rate coefficients Tj, which are the rates for 
dipole orientations parallel to the principal axes Vj, lead- 
ing to: 

r(e d ) = ^(e d )r min + /? 2 2 (e d )r mcd + (3 2 3 (e d )T m ^. (5) 

Equation ([5]) describes the emission-rate surface as a 
function of dipole orientation r(e<i), which is a central 
result of our work. The emission rate coefficients Ti are 
equal to: 

r* = -r 2vf.ImG. Vi =- l9ll (6) 

and are via G functions of the frequency and the dipoles' 
position: Ti = Ti(u>,r). Assuming known principal rates 
Ti, the emission-rate surface r(ed) is always a quadratic 
form on the unit sphere. Moreover only quadratic forms 
of signature s = J2(sign(Ti j) = 3 can occur 20], since 
emission rates are physically constrained to be positive 
for all orientations. Therefore, polar plots of the rate 
versus dipole orientation - henceforth called emission- 
rate surface - take on only specific shapes classified by 
the ratios of r m i n ,r mo d, r max , with three perpendicular 
symmetry axes, regardless of the nanophotonic system. 
We remark that while Eq. J5| may appear as the defining 
equation of an ellipsoid, the emission-rate surface is not 
an ellipsoid since the problem is not about calculating 
a level surface of Eq. ([5]), which would be equivalent 
to constraining to yield a fixed T in Eq. ([5]), rather 
than constraining to the unit sphere. Our result that 
emission surfaces are always necessarily quadratic forms 
defies the intuition (as sketched in Fig. [T|) that emission 
rates inherit the symmetry of the nanophotonic system. 

Regarding the assumptions we require to arrive at the 
quadratic form for the emission rate surfaces, we note 
that we have assumed real dipole moment in Eq. (fol- 
lowing Ref. [J]) and that we used reciprocity to ensure 
real and symmetric Im(G(r, r)). In case of reciprocal 
media it is easy to show that our results are also valid 
for complex transition dipole moments, and not just for 
real dipole moments. Furthermore, if we assume a real 
dipole moment, it appears that our results are also valid 
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FIG. 2: (color) (A) Drawing of a two-level quantum emitter at distance h above a mirror. (B) Emission rate versus scaled 
distance (wave vector times distance kh) for a dipole perpendicular and parallel to a perfect mirror ||. (C) Three-dimensional 
surfaces representing the orientation dependent spontaneous emission rate in real space, (i) One maximal emission and two 
equal minimal rates give a peanut-shape (at kh = 0.4 in (B)). (ii) Two equal maximal rates and one minimal rate give an oblate 
spheroid (kh — 2.3 in (B)). (iii) Three equal maximal rates give a sphere (kh = 3.2 in (B)). (D) Most general shape when all 
principal rates are different (r max > r mo( j > r m i n ) and the principal axes are rotated from the (x, y, 2)-axes. Color scales are 
linear from r m j, n to F max (colorbar in (C)). 



for metamaterials that violate reciprocity, i.e., in case 
Im(G(r,r)) is not symmetric or even not diagonalizablc. 
Since Im(G(r, r)) is still real it will nonetheless give rise 
to a quadratic form that can be transformed to a prin- 
cipal axis system [20| . The physical requirement that 
rates are positive for all dipole orientations furthermore 
ensures that the signature of the quadratic form remains 
3 even in the nonreciprocal case. 



B. Generic shapes of the emission-rate surface 

Figure HKC,D) categorizes all possible shapes of the 
emission rate polar plot. Fig. djC) is relevant for the 
mirror, with principal axes parallel (x, y, degenerate) and 



perpendicular (z) to the interface. Fig. HjC(i)) shows 
the emission-rate surface for the case where emission 
is enhanced along a single dipole orientation r max 3> 
Tmin = Timed- This situation appears at a reduced dis- 
tance kh = 0.4 close to the mirror. Here, the emission- 
rate surface looks like a highly anisotropic peanut, con- 
stricted to a radius r m i n in the x, y-plane, and extending 
to r max along the z-axis. Fig. [5JC(ii)) shows the ori- 
entation dependent emission rate for a single inhibited 
axis with r m i n < r me( j = r max , at kh = 2.3 near a mir- 
ror. Qualitatively, the emission-rate surface resembles 
an oblate spheroid; when the minimum rate is much less 
than the other two rates (see Fig. below), the surface 
develops a concave indentation with a donut-like shape. 
Fig- HlC(iii)) shows the emission-rate surface when the 
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rate is equal along all three main axes (kh = 3.2). The 
emission-rate surface is simply a sphere, as it is in any 
isotropic homogenous medium. 

Fig. \T[D) shows the emission-rate surface for the most 
general case when i) the rates along the main axes are 
all different (r m ; n < r mo d < r max ), and ii) the princi- 
pal axes v m i n)ln ed,max have an arbitrary orientation with 
respect to the laboratory frame. Clearly, the emission 
rate is not extremal for a dipole parallel to any of the 
(x, y, z)-axes. An important feature of the emission-rate 
surfaces is that they allow for an easy inspection of both 
the anisotropy of the emission rates, and of the favorable 
dipole orientations compared to the usual (a;, y, z)-axes 
in real space. 



III. EFFICIENT METHOD TO CALCULATE 
EMISSION-RATE SURFACES 

In many cases of practical interest, neither the Green's 
function G nor the principal axes {vi} are a-priori 
known. Often algorithms based on a summation over all 
photon modes are used that only yield the rate T for tar- 
get orientations chosen as a priori input. Reconstruct- 
ing emission-rate surfaces as in Fig.[2]by a dense sampling 
of orientations is not viable with such algorithms, due 
to prohibitive computation times. A poignant example 
is the calculation of emission rates in photonic crystals 
that requires a summation over up to 10 6 Bloch modes, 
the calculation of each of which requires diagonalization 
of a 10 3 x 10 3 matrix, even for a single dipole orienta- 
tion [ll|, [2l[ . A popular alternative method that can con- 
veniently yield the emission rate for a single orientation 
is the finite difference time domain (FDTD) simulation 
method (22J. However, it appears difficult to calculate 
off-diagonal elements of the Green tensor. Since the var- 
ious field components are not calculated on identical grid 
points, FDTD does not truly yield a Green dyadic on a 
well-defined position r. Hence, even if an algorithm is 
known to calculate rates at fixed orientations, it is un- 
clear how to find the principal axes and rates, since Im(G) 
is simply not available for diagonalization. In view of the 
computational cost of evaluating the radiative rate at a 
single dipole orientation, the main problem is to find out 
for how many and for which orientations the emission 
rate must be calculated to completely and exactly char- 
acterize the emission-rate surfaces. Here we describe an 
efficient method to find principal emission rates and ori- 
entations by evaluating the LDOS at the least possible 
number of input orientations. 

We use the well-known fact that any function on the 
unit sphere is conveniently expanded in spherical har- 
monics Y lm (9, <f>) — P/ m (cos(#))e"™*. Since the emis- 
sion surface is a quadratic form, we can apply the well- 
known fact that all quadratic forms on the unit sphere 
can be represented exactly by an expansion containing 



only terms up to I — 2, so that 
2 / 

r(eJ = ^^ flI X(cos(«))e'^ (7) 

i=0 m=— I 

An easy proof that no terms beyond I = 2 are needed is 
obtained by expressing the spherical harmonics in terms 
of cartesian coordinates, rather than polar coordinates 
on the unit sphere 20], or conversely by expressing the 
coefficients /3i in terms of polar coordinates relative to the 
{vj} axis system. This substition leads to a trigonomet- 
ric expansion for r(e,j) with terms that are quadratic in 
cosines and sines of 9 and of </>, see Appendix IA1 Eq. (jAip . 

The expansion coefficients for the spherical harmonic 
expansion are given by inner products 

a lm = (T{e d ),Y lm ) = / #/ d6r(e d )Y lm (6,<t>)sm(6), 

JQ JO 

(8) 

similar to the coefficients appearing in discrete Fourier 
transformations, but now for transformation on the 
unit sphere. Mohlenkamp has developed a fast Fourier 
transform method to calculate the coefficients numeri- 
cally [23( | . which requires a sampling of rates T at a dis- 
crete set of orientations, similar to the numerical evalu- 
ation of discrete Fourier coefficients by the sampling of 
a periodic function on a discrete set of points. In this 
approach, the integral expression ([5]) for the expansion 
coefficients for expanding a function / is replaced by a 
discrete weighted sum: 

ai m = 2J Wkf(&k, <t>k)Yi m {9k, <t>k) sin(8 k ), (9) 

k 

where k runs over the finite set of sampling points. Such a 
discrete approximation to the expansion coefficients a\ m 
is in fact exact for all functions / that are exactly equal 
to a finite series of spherical harmonics up to order Z max 
if: i) the angles (9k,4>k) are chosen as the roots of the 
basis functions of order I = Z max + 1, and ii) the Wk 
are appropriate weights. In the present case ^ max = 2. 
Thus the special points are the 18 roots of the spherical 
harmonics Y^ m (9,(j)) of order I = 3. Furthermore, one 
may appreciate that the spherical harmonic transform is 
a simple Fourier transform over <^>, and a Legendre trans- 
form over cos 9. The weights Wk are hence the weights 
appropriate for Gauss-Legendre quadratures of order 3. 
Explicitly, the 18 special points occur at azimuthal an- 
gles <fi = mir/3 (m = 0,1,... 5) and at polar angles 
9 = arccos(-\/3/5), 7r/2, arccos(— ^/3/5). The weights w 
only depend on 9, and are 5/9 for 9 = arccos(±-\/3/5) 
and 8/9 for 9 = tt/2. Since one half of the 18 points (see 
figure [3]) is antipodal to the other half, inversion symme- 
try of the emission rate means that the rate need only 
be evaluated for 9 dipole orientations in order to find the 
full spherical harmonic expansion. 
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FIG. 3: (color online) Special orientations, i.e., points on the 
unit sphere, for which the decay rate needs to be calculated 
in order to fully reconstruct emission-rate surfaces. The blue 
and red points together are the roots of I = 3 spherical har- 
monics, corresponding to (f> = mn/3 (m = 0, 1, 2, ... 5) and 
8 = arccos(\/3/5), 7r/2, — arccos(\/3/5). Due to inversion 
symmetry, rates are equal for antipodal orientations. This 
makes calculations for half the points (e.g., the blue ones) su- 
perfluous, leaving 9 distinct orientations (red) for which rates 
must be calculated in order to find principal rates and orien- 
tations. Note that these points do not have equal weights in 
Eq. (O (weight is 8/9 for points on the equator and 5/9 for 
other points). 



IV. RESULTS FOR SEVERAL 
NANOPHOTONIC EXAMPLES 

To illustrate our analysis, we discuss the emission dy- 
namics of a quantum emitter inside a photonic crystal, 
illustrated in Fig. [4jA,B). These complex systems have 
extreme variations of the emission rate versus frequency 
on account of a bandgap where emission is completely 
inhibited @. To obtain the rate for an emitter of arbi- 
trary orientation in a Si inverse opal, we have calculated 
the LDOS for the 9 special orientations by summing over 
all Bloch eigenmodes [24j . The crystal has a first order 
'pseudogap' at reduced frequency 0.55, and a photonic 
bandgap from 0.852 to 0.891. Figure [4jB) shows the 
emission rate for a salient position in the unit cell (c/. 
Fig. BJA)): the rate is anisotropic for frequencies near 
the pseudogap, since it differs for dipoles pointing in ci- 
ther x, y or the z direction, which are the cubic symmetry 
axes of the crystal. One might be tempted to perceive 
the behavior to be as simple as a mirror, since it is the 
same for both x and y. However, a plot of the maxi- 
mum, medium, and minimum emission rates (Fig. HJC)) 
shows that this perception is completely wrong: Already 
at low frequency up to the pseudogap, the emission rate 
is strongly anisotropic. While anisotropic behavior in 
the long-wavelength limit may seem surprising, its origin 
in electrostatic depolarization effects has been discussed 
before [26|, [27j ■ The maximum rate occurs for dipole ori- 
entation e c i = {— 1, 1,0)/ v2, and is much larger than the 
rate for any of the x, y, z orientations, whereas the mini- 



mum rate for = (1, 1, 0)/V2 is much smaller. At high 
frequency (a/A > 0.6) up to the bandgap, the orienta- 
tion of maximum rate changes to = (0,0, 1). While 
it is clear from Fig. 0JG) that the orientation dependent 
emission rate is much more complex than expected from 
(B), Fig. \4^C) hardly gives an intuitive picture of the 
orientation-dependent behavior. 

Therefore, we plot in Figure [4^D) emission-rate sur- 
faces versus frequency. At frequencies below the pseu- 
dogap, the emission-rate surface is peanut-like, reveal- 
ing that the emission rate is high for a "horizontal" 
dipole orientation, and inhibited for the 2 perpendicu- 
lar orientations. At the pseudogap, the emission-rate 
surface suddenly changes to donut-like, since the rate is 
high for two orientations and low for a third orienta- 
tion. At even higher frequencies, the emission-rate sur- 
face becomes again peanut-like - with donut-like behav- 
ior near 0.8 - but with a different orientation than be- 
low the pseudogap. The maximum emission rate is up 
to 20-fold enhanced, and the anisotropy (r max /r min ) is 
strong with peaks up to 340. In this particular exam- 
ple, the high symmetry at this spatial position fixes all 
principal axes. To demonstrate the applicability of our 
method to general, nonsymmetric, cases we have also 
studied low-symmetry positions at constant frequency, 
see Fig. UA). Again strong anisotropies occur, with the 
maximum-emission axis (or inhibition-axis) continuously 
changing direction as a function of source position. We 
conclude that emission-rate surfaces provide a compact 
representation of the rich behavior of the dependence of 
the emission rates on dipole orientation. 

We emphasize that our classification of emission dy- 
namics by means of emission-rate surfaces is by no means 
restricted to dielectric systems and can also be applied 
to dissipative nanophotonics systems that are of modern 
interest, such as plasmonic and metamaterial structures. 
Our analysis rests purely on the symmetry of the Green 
dyadic in Eq. which in the presence of optical ab- 
sorption describes the total decay rate (radiative rate 
plus induced nonradiative rate) of a quantum emitter. 
As an example, we discuss the textbook case of an emit- 
ter near a plasmonic sphere [32l . HH, using the known 
Green's function [H (c/. Fig.^A)). Figure^C) shows 
that the emission-rate surface for the total decay rate has 
a donut-like shape (r m i n < r mcc i = r max ) with 16-fold 
enhanced rates for a dipole parallel to the surface, and 
5-fold enhanced for a perpendicular dipole. For a fixed 
dipole orientation 0, EH , the angular distribution of the 
radiated power reveals a well-known five-lobed structure 
(B). A comparison of (B) and (C) illustrates the main 
differences between radiation patterns and emission-rate 
surfaces: radiation patterns are relevant to a single dipole 
orientation and do not necessarily have any symmetry, 
or are free to follow any symmetry inherent in the envi- 
ronment. Emission-rate surfaces on the other hand are 
relevant to all orientations and have a symmetry limited 
by the quadratic form. 
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a/A a/A 

FIG. 4: (color) Emission rate for a quantum emitter in a photonic bandgap crystal. (A) Left: l/8th of a cubic unit cell, blue 
dashed lines delimit the primitive cell. Right: emission-rate surfaces on x = 0.2, y = 0.3, variable z (red line in left panel) 
at reduced frequency a/A = 0.94 (a is lattice parameter). Surfaces are colored by relative rate (scalebar on right), and have 
constant size. (B) Emission rate for dipole centred in a window of a Si inverse opal [r = 1/4(1, 1,0), black dot in (A)] with 
orientations = (1,0,0), (0,1,0) (black curve), and (0,0,1) (blue dashed-dotted curve) versus a/A. The rate is normalized 
to the one in vacuum. The blue vertical bar indicates the photonic bandgap. Inset: cubic unit cell. (C) Maximum, medium, 
and minimum emission rates T m ax (+), r mec j (x), Y m in (|) compared to rates for orientations = ( — 1, l,0)/v2 (full curve), 
(0,0,1) (dashed-dotted), (1, l,0)/v2 (short dots). (D) Emission-rate surfaces at select frequencies show strong changes in 
shape. The size of the surfaces is in proportion to the absolute emission rates, and colorscales range from Y m i a to r max . Dashed 
curve: F max . 




FIG. 5: (color) (A) Drawing of a two-level quantum emitter at 20 nm distance from a plasmonic Ag sphere [23] with radius 
R = 80 nm. (B) Angular distribution of the radiated power versus solid angle Q for a single dipole orientation parallel to 
the surface; the pattern has a complex five-lobed structure. (C) Emission-rate surface showing the emission rate versus dipole 
orientation. The pattern has a donut-like shape. 
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V. DISCUSSION 

Since the analysis in this paper is based on Im(G(r, r)) 
it is strictly valid for the total decay rate modification in- 
duced by the nanophotonic environment. Explicitly, in 
the case of losses our proof only holds for the sum of the 
radiative rate and the non-radiative rate (r ra d + r nonra d), 
and not for the radiative rate r ra( j separately. To analyze 
the radiative emission-rate surfaces one would need to an- 
alyze the far-field integral of the radiated power (quantity 
in Fig. GD(B)) as a function of the source orientation. A 
priori it is not at all clear that such radiative rate sur- 
faces need have a quadratic form. Indeed, we have not 
succeeded in proving the quadratic form for the radiative 
rate in the lossless case by analysis of far-field integrals, 
i.e., without identifying r ra d = r to t and subsequently 
analyzing Im(G(r, r)). We have numerically calculated 
radiative emission rate surfaces for many low-symmetry 
dissipative plasmon sphere clusters, and have not found 
any example in which the radiative emission rate surface 
was not quadratic. Although a rigorous proof is beyond 
the scope of this paper, we therefore anticipate that the 
quadratic form not only holds for total decay rates, but 
also for radiative decay rates. 

A class of quantum emitters with a single transition 
dipolc moment arc fluorescent molecules, such as laser 
dyes Q. For such emitters, emission-rate surfaces can be 
observed if their orientation is controlled, e.g., by attach- 
ing them to liquid crystal molecules that are oriented in 
external fields 31]. If one can tune the orientation of 
an emitter, this opens a novel opportunity to "switch" 
spontaneous emission from inhibited to enhanced and 
vice versa. The emission-rate surfaces reveal that optimal 
switching always requires a dipole rotation by 90°, since 
minimal and maximal emission rates always occur along 
the mutually perpendicular main axes. Alternatively, one 
could tune semiconductor nanowires with oriented dipolc 
moments. For self-assembled and colloidal quantum dots 
with dipoles in a x' , y' plane, we expect to probe the x 1 , y' 
cross-sectional average of the emission-rate surface of the 
relevant nanophotonic system. 

Since arbitrary orientations do not usually coincide 
with principal dipole orientations, most prior work on 
specific systems has been incomplete, since no principal 
rates has been reported. While such incompleteness does 
not affect the orientation averaged rate (see Appendix [A"|) . 
it does affect the understanding of dynamics of orienta- 
tional dipole ensembles @, [T3, EH • Such a decay is a sum 
of single exponentials with a rate distribution given by 
the emission-rate surface. Any observable derived from 
time-resolved decay beyond the orientation-averaged rate 
(Tr(Im(G))) requires knowledge of the principal rates, 
which is thus relevant to many physical situations in 
nanophotonics . 

In classical optics, the imaginary part of the Green 
dyadic is not only relevant for radiating dipoles. Indeed, 
the imaginary part of the Green dyadic has also been con- 
nected to the so-called coherency matrix (or the electric 



cross-spectral density tensor) (34j for black body radia- 
tion. In general, the 3x3 coherency matrix describes 
second-order spatial correlations of the electric field, and 
can be understood as a generalization of Stokes param- 
eters to quantify the polarization of near fields locally. 
Within this framework, a description of local polarization 
by polarization ellipsoids directly points at a quadratic 
form of the coherency matrix, since ellipsoids are level 
sets (rather than polar plots) of an equation of the form 
in Eq. ([5]). It should be noted that the coherency matrix 
depends on the incident source that generates the local 
electric field. In the particular case that the field is due 
to black body radiation the coherency matrix reduces to 
the imaginary part of the Green dyadic Im(G(r, r)), as 
derived by Setala et al. [35| . However, it is important 
to realize that for this identification of Im(G(r,r)) with 
the coherency matrix to hold, the source is required to 
be a statistically homogeneous and isotropic distribution 
of radiating currents, and the medium is supposed to be 
non-dissipative [H| . This is diametrically opposite to the 
analysis of spontaneous emission sources presented here, 
which concerns localized and oriented sources and is valid 
without limitation on material dissipation. It is exciting 
that our method to find principal rates and orientations 
can be directly adapted to calculate the local polarization 
properties of black body radiation. 



VI. SUMMARY 

We have theoretically studied the spontaneous emis- 
sion rate of a two-level quantum emitter in any nanopho- 
tonic system. We derive a general representation of 
the dependence of emission rates on the orientation of 
the transition dipole by only invoking symmetry of the 
Green function. The rate depends quadratically on ori- 
entation and is determined by rates along three princi- 
pal axes. We show that these principal rates and axes 
can be easily calculated without evaluation of the full 
Green function. Furthermore we show that visualization 
of emission-rate surfaces as determined from principal 
rates provides great insight on how preferred orientations 
for enhancement (or inhibition) depend on emission fre- 
quency and location, and on strategies to actively switch 
emission rates by the dipole orientation, as shown for a 
mirror, a plasmonic sphere, or a photonic bandgap crys- 
tal. 
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APPENDIX A: DISCUSSION OF AVERAGE 
EMISSION RATE 

A remarkable fact is that the orientation-average rate 
(F) can always be calculated from the LDOS at just three 
perpendicular orientations, which need not coincide with 
the principal axes {v m i n , v mc d, v max }. First, we calcu- 
late the orientation averaged rate by integration over the 
full emission surface. Without loss of generality we align 
x,y,z with the principal axes, so that the orientation- 
dependent rate is: 



Integration over the full emission surface clearly shows 
that the orientation-averaged emission rate is equal to 
the mean of the three principal rates, and hence (r) = 
(Trd 2 u>/heo)(2uj/irc 2 ■ Tr(Im(G(r, r, u>))). The invariance 
of the trace of any matrix under arbitrary basis rotation 
implies that the average rate in Eq. (|A2p can be calcu- 
lated from the rates at any randomly chosen but mutually 
orthogonal directions x, y, z as 



T(9, (f>) = r min cos 2 4> sin 2 6>+r mod sin 2 sin 2 6»+r max cos 2 6. 

(Al) 

By straightforward integration, the orientation-averaged 
rate (r) is 



(r} = -(r 2 , + r. 



(A3) 



(r) 



i 

4^r J 



i sin 6d6 



1 



(r 



minTi mcd 



(A2) 
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